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Abstract 

An on-shell renormalisation programme for the chargino/neutralino and the sfermion 
sectors within the Minimal Supersymmetric Standard Model as implemented in a 
fully automated code, SloopS, for the calculation of one-loop processes at the collid- 
ers and in astrophysics, is presented. This is a sequel to our study in Ref. [T] where 
an on-shell renormalisation of the Higgs (and the gauge/fermion) sector is performed. 
The issue of mixing is treated in a unified and coherent manner in all these sectors, 
in particular we give some new insight into the renormalisation of the mixing angle 
in the sfermion sector and like with the Higgs sector and the issue of tan j3 we discuss 
different schemes. We also perform numerical comparisons between our code SloopS 
and different results found in the literature. In particular we consider loop correc- 
tions to the neutralino and sfermion masses, chargino pair production and stau pair 
production in e + e~ colliders, as well as a few decays of the heavier chargino. For all 
these observables, we analyse the tan (5 scheme dependence using different definitions 
of this parameter and comment on the impact of using different renormalisation of 
the mixing parameter in the sfermion sector. 



1 Introduction 



The description of the Higgs within the Standard Model is unsatisfactory as it poses the 
problem of naturalness. Besides, the Higgs particle is still missing. Moreover there is 
overwhelming evidence that there is a large amount of Dark Matter that can not be ac- 
counted for by any of the particles of the Standard Model, SM. All this points to New 
Physics. The best motivated model of this New Physics is undoubtedly supersymmetry 
that rests on solid theoretical grounds and allows for full calculability and therefore predic- 
tions. Full calculability is not, by itself, a sacrosanct virtue but it must be admitted that 
supersymmetry addresses some of the problems of the SM. Indeed, although the primary 
motivation for supersymmetry as implemented in the MSSM, Minimal Supersymmetric 
Standard Model, was to solve the hierarchy and naturalness problem it was soon realised 
that the model contained an excellent candidate for cold dark matter besides incorporating 
almost naturally the gauge unification. However, predictions of the MSSM based on tree- 
level calculations predict a Higgs that is lighter than the Z mass. By now this is ruled out. 
It is only through radiative corrections that the MSSM has survived. Radiative corrections 
are therefore essential. Moreover the next generation of experiments at the colliders will 
reach unsurpassed precision which will need computations beyond the tree approximation. 
Extraction of the cosmological parameters that are used to measure the relic density of 
cold dark matter have recently reached an accuracy that will also soon compete with the 
accuracy we have been accustomed to from the LEP era. Precision loop calculations within 
the MSSM are therefore a must. It must be said that quite a lot of these calculations have 
been performed, even though the bulk of these have been made for collider observables 
and indirect precision measurements such [g — 2)^,6 — > 57,... Very little has been done 
concerning the cross sections relevant for dark matter annihilation that enter, for example, 
a precise prediction of the relic density. It rests that these calculations have been done 
piecemeal and quite often within different renormalisation schemes. 

One of the reasons that these calculations have been done piecemeal is that the MSSM, 
though minimal, still contains a large number of particles and a very large number of 
parameters especially through the soft-susy breaking terms for example. This explains 
why different groups have concentrated on different sectors of the model. Performing loop 
calculations with so large a number of parameters and huge number of interactions is an 
almost untractable task especially if one has to be ready to perform precision predictions 
for any process or at least a large number of processes as it occurs for example with the 
calculation of the relic density where many processes and sub-processes are at play for 
a particular choice of parameters. One has to rely on a fully automatised code for such 
calculations. SloopS is such a code with an automatisation starting already from the 
implementation of the model file. Instead of coding by hand all the Feynman rules which 
usually constitute the model file and realising that for one-loop applications one needs to 
also enter the full set of counterterms, SloopS relies on a much improved version of LanHEP 
[2] to automatically generate the model file. Through LanHEP one writes the Lagrangian 
in a compact form through multiplets and the use of the superpotential. The improved 
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version of LanHEP has built-in rules for shifting fields and parameters thus easily generating 
the set of counterterms. This approach therefore takes care of generating the few thousand 
Feynman rules for all the vertices needed for the calculations of any one-loop process in 
the MSSM. 

The model file thus generated is interfaced to the bundle of packages FeynArts [3], FormCalc 
[I] and LoopTools [5], that we will refer to as FFL for short. This code has recently been 
used very successfully for the first calculation of a number of processes that enter the pre- 
diction of the relic density of Dark Matter [6] as well as some one-loop induced processes 
of relevance for indirect detection [7]. 

The aim of the present paper is to first give some details on the renormalisation scheme 
that is implemented in SloopS and in particular how the sfermion sector and the neu- 
tralino/chargino sector are treated. This is a follow up to our paper detailing the renor- 
malisation of the Higgs sector where apart from the implementation of the scheme we 
brought up crucial issues related to the definition of tan (3, the issue of gauge invariance 
and the impact of different schemes on observables in the Higgs sector. The present pa- 
per will also compare one-loop predictions in the sfermion and chargino/neutralino sector 
based on different schemes for tan (3. We will also make some interesting observations and 
analyses concerning the treatment of mixing in these sectors, especially how one could 
define a process independent mixing angle in the sfermion sector. 

The paper is structured as follows. In Section [21 we give a brief summary of the renormal- 
isation scheme used in the code for the Higgs sector and the SM-like sector that includes 
the gauge and fermion parts. In the same section we also present a general overview 
of our approach. Section [3] deals with the sfermion sector, both squarks and sleptons, 
that we use in SloopS. In Section H] we detail our on-shell renormalisation scheme in the 
chargino/neutralino sector and comment on some alternatives for the choice of the in- 
put parameters. Section [5] illustrates the use of the code for some applications. We will 
give results for the one-loop corrections to the masses of the heavier neutralinos and the 
sfermions that are not used as input in our schemes. We also present results for the 
one-loop calculation of chargino pair production and sfermion pair production at a linear 
collider, e + e~ — ► XiXi an d e + e~ — > fjfj comparing. Finally we compare our results with 
those of Grace-SUSY [8] taking as examples a few decay channels of the heavier chargino 
for a certain choice of parameters. In all these examples the tan/3-scheme dependence 
is also studied thus complementing the scheme dependence that we studied for observ- 
ables within the Higgs sector and for annihilation processes of interest for the relic density 
computations. Section [6] gives a brief summary and outlook. 
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2 Renormalisation: The general approach, the gauge, 
the fermion and the Higgs sector 



Our renormalisation of the MSSM, with CP conservation with all parameters taken real, 
follows the same strategy and the same procedure that we adopted for the renormalisation 
of the Standard Model, see [9]. In particular we strive for an on-shell renormalisation of 
the physical parameters. Counterterms to these parameters are gauge independent. Wave 
function renormalisation is introduced is order that the residue of the two-point function, 
the propagator, is unity for the physical state on its mass shell, as well as to eliminate 
any mixing between the physical fields when these are on-shell so that the qualification 
as a physical field is maintained order by order. Naturally, these field renormalisation 
constants are not needed if one only requires that the observables of the S'-matrix are 
finite but one does not insist that all the Green's function to be finite, see [9]. On the 
technical side this field renormalisation avoids that one includes in the calculation of matrix 
elements loop corrections on the external legs. Moreover there is no need to consider field 
renormalisation for the unphysical fields like the Goldstones bosons or on the current fields 
before mixing. Talking about the Goldstone fields a very powerful feature of SloopS is 
the use and implementation of a non-linear gauge fixing condition [TQl [1] . The gauge- 
fixing condition furnishes eight gauge parameters (a, ft, 5,£u, R, p, e, 7) on which we could 
perform gauge parameter independence checks, beside the ultraviolet finiteness checks. 
The gauge-fixing writes 

C GF = -—F + F--—\F z \ 2 -—\F^\ 2 ,with 
iw 2£z 2£ 7 

F+ = (d^-iea 1 ^-ie—pZ^W^ + + i£ w — (v + 5h + uH + tpA + ikG )G + , 

Sw 2s\y 

F z = d »zl + (i z —{v + ~eh + ; 1 H )G\ 

F 1 = drf. (2.1) 

As extensively stressed in [9] and [1] the gauge fixing term is considered renormalised. h° 
and H° are, respectively, the lightest and heaviest CP-even Higgses, A is the CP-even 
Higgs, G 0,± are the Goldstone bosons and, W ± , Z, 7 are, with obvious notations, the gauge 
fields. We have c\y = cos 6\y = M w /M^j . We work with £w,z,i = 1 in order not to have 
to deal with too high a rank tensors concerning the loop libraries, see p]. 
Another crucial feature of our renormalisation program is our treatment of the mixing 
which occurs in all sectors of the MSSM. In general, fields are expressed in the current basis. 
They, however, mix. Physical mass eigenstates fields are obtained from these current fields 
through some rotation matrix at tree-level. We consistently take, in all sectors, this matrix 
to be renormalised and therefore no extra counterterm is introduced to this matrix. At 
one-loop, this will still leave some transitions between fields, however field renormalisation 

*To avoid clutter we use some abbreviations for the trigonometric functions. For example for an angle 
9, cos 9 will be abbreviated as eg, etc... tp will then stand for tan/3. 
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is defined to precisely get rid of any residual mixing when the physical particles are on-shell. 
Therefore inducing counterterms for the rotation matrix is redundant and not helpful. 
Let us now briefly recap on the renormalisation of the gauge, fermion and Higgs sector. 

2.1 The fermion and gauge sector 

The fermion sector as well as the gauge sector are renormalised on-shell. It means, for 
example, that the gauge boson masses My/± and M^o are defined from the pole mass, 
imposing the one-loop on-shell condition on the mass counterterms as 

5M^ ± = -Re^ ±w± (M^ ± ) , 5M 2 Z0 = -#e£|o z o(Mf ) . (2.2) 

The electric charge e is defined in the Thomson limit. Since MSSM processes and pa- 
rameters are taking place at the weak scale, the effective gauge coupling constant is of 
order a(M| ) which includes large logarithms from the very light standard model charged 
fermion masses. It is useful to reparameterise the one-loop corrections in terms of this 
effective coupling in order to absorbs these large logarithms as we will see later. 

2.2 The Higgs sector 

The renormalisation scheme and renormalisation procedure at one-loop in the Higgs sector 
that we adopt in the code is detailed in Ref. p]. The only ingredient that makes its way 
from the Higgs sector and the Higgs observables to the chargino/neutralino sector and the 
sfermion sector is the ubiquitous tp and its renormalisation. We use the same notation 
as in [9]. At tree-level, tp is defined by the ratio of the two vacuum expectation values 
tp — V2/V1. At one-loop, as pointed out in Ref. [TTl [1] it is difficult to find a proper 
definition for tp. In [I] we critically discussed the issue of gauge invariance as regards 
different definition of tp and looked quantitatively at the scheme dependence introduced 
by tp in some Higgs observables. We will extend this investigation in our applications to 
observables involving the sfermions and the chargino/neutralinos. We therefore consider 4 
definitions which are detailed in [I]. 

• v4 rT -scheme. 

tp is extracted from the decay A — > t + t~ to which the QED corrections have been 
subtracted. This leads to a gauge-independent counterterm. In Ref. [12] the decay of 
the charged Higgs boson H + into r + and associated neutrino v T has been suggested. 
This would qualify as a gauge independent definition, the advantage of our A 
is that the full QED corrections can be extracted most unambiguously. 

• Mif-scheme. 

Here the heaviest CP-even Higgs mass M H o is taken as input. This definition is 
obviously gauge independent and process independent but unfortunately, we remarked 
that it induces large corrections in many cases. 
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DR-scheme. 

Here only the ultra-violet part of an observable such as A — > t + t~ (or any other 
definition but within the linear gauge, see p]), is extracted. In this scheme, the 
counterterm depends explicitly on a scale JL. This scale Jt is fixed at M40. 

DC PR-scheme pj]. 

dtp is extracted from the A°-Z° transition at q 2 = M A0 , 

r. DCPR 

-r 1 = - T , ReZ A o z0 (M 2 AO ) . (2.3) 

The self-energy of the A — Z° transition at large t/3 is dominated by the bottom/tau 
loops because of the A°bb vertex which is proportional to m^tp and thus enhanced 
when t/3 becomes large, 

j£ dcpr t a 2 1 

~T ^-^^^^(3mtB (M%,mlm 2 b )+m 2 T B (Ml ,ml,m 2 T )) . (2.4) 

tp S 2 (3 C W M Z 47T 

The loop functions Bq is defined in [14] . At large tp S2/3 ~ the finite part of Stp/tp 
in the DCPR scheme is of order ti. This scheme is not gauge independent and would 
depend on some parameter of the non-linear gauge for example. When comparing the 
results of observables within this scheme we will set all non-linear gauge parameters 
to zero, i.e. we will be specialising to the linear gauge. 



3 The sfermion sector and its renormalisation 

The sfermion sector comprises the superpartners of the fermions of the Standard Model 
where the interaction fields are the chiral left and right states. We do not consider gener- 
ation mixing. For each generation, the field content is therefore the doublet Ql = (ul, g?l) 
and singlets ur and d,R for the squarks. For the sleptons we have E L = (u L , e L ) and e^. In 
case the corresponding Yukawa coupling is zero with vanishing fermion masses, we expect 
no u L — u R and — d R mixing, so that the physical fields are Ul.Ur and dL,d R in the 
squark sector. Let us briefly recall where the mass parameters of the sfermion sector orig- 
inate from, and how many can be identified and defined solely within the sfermion sector, 
once for example the Higgs sector and gauge sector have been identified. 

• The soft supersymmetry breaking terms 

C lft = z2 ], li<f> fi = QL,L L ,u R J R ,~e R (3.1) 

fi 



e l3 (y^A u HiQiu R + y^A d H{Qid R 



+ ^l AeH f L {e R + h.c.). 
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Our conventions for the Higgs doublet and the vacuum expectation values of these are 
defined in pQ. Supersymmetry breaking therefore provides the soft scalar masses M~ 
Eq. (13.11) and the tri-linear scalar coupling Af parameters Eq. (13. 2p . for / = e, u, d of 
one generation. The contribution of the latter vanishes in the chiral limit where the 
mass of the fermion, rrif, vanishes. The latter generates not only a contribution to the 
mass of the different sfermions but also contributes to the coupling of the sfermions 
to Higgses and Goldstones. As known, because of the SU(2) symmetry, there is only 
one soft mass parameter for the up and down left component of the scalars. 

• Sfermion masses get also a contribution from the usual Yukawa mass terms, these are 
proportional to the corresponding rrif. 

• We also get contributions from the supersymmetry conserving F terms. The F(f) 
contribution does not mix left and right explicitly (though it is proportional to the 
corresponding fermion masses, rrij). This only generates couplings to Higgses. The 
F(Hi 2) involve the /x parameter and generate supersymmetry conserving tri-linear 
scalar coupling. They lead to left-right mixing which is proportional to ruffj,. 

• There are also D term contributions, chirality conserving, proportional to the gauge 
boson masses. These give contributions to the sfermion mass terms, //, Higgs cou- 
plings ffH, G and quartic scalar couplings: //// and ffHH. Once the gauge and 
Higgs sector have been renormalised these contributions are also. 

These simple observations show that since the Af terms and ji contributions do not act simi- 
larly on the mass term and the Higgs couplings of sfermions, renormalisation of the sfermion 
two-point functions (mass, mixing and wave function renormalisation) is not enough to 
completely renormalise processes with ordinary standard particles and the sfermions. One 
needs also to define a renormalisation to the \x parameter. This is most conveniently done 
from the chargino/neutralino sector. Note however that the Higgs coupling to sfermions, 
can provide an alternative definition to fi. 

3.1 Renormalisation of the Squark sector 

We show in detail the different steps specialising to those squarks with mixing. The case 
with no-mixing is then trivial. 

3.1.1 Fields and parameters at tree-level 

The tree level kinetic and mass term for the squarks q = u,d are given by, 

& = ~\ ( Ml « ) ( £| ) + U ?R ) M\ ( | ) , (3.3) 
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with the 2x2 non- diagonal mass matrix 



M 2 qLL M? LR 
M? LR M? RR 



(3.4) 



The different components of this matrix are, 



M 2 LL = Ml L + m 2 q + c w (Tl-Q q s 2 w )M 2 z , (3.5) 

Ml RR = Ml + m 2 q + c 2 pQ q s 2 w M 2 Zj (3.6) 

-2T„ 3 N 

L qLR — " L q\^q ~ V« j3 



Ml R = m q {A q -iit B (3.7) 



M% is the soft-SUSY-breaking mass parameter of the 577(2) l doublet, whereas M| is the 
soft-SUSY-breaking mass parameter of the singlet. T q and Q q are the third component 
of the isospin and the electric charge respectively. M~ LR is the mixing parameter that 
has contributions from both the higgsino supersymmetry conserving mass parameters and 
the tri-linear supersymmetry breaking term. This induces mixing between the left and 
right components. This mixing vanishes for sfermions associated to massless quarks but is 
important especially for the third family squarks. Note that this mixing can also vanish, 
at tree-level, even for massive quarks for exceptional A t = ix/tp for stops and = fitp for 
sbottoms. 

If fi is to be determined from the chargino/neutralino sector, this sector involves 5 new 
parameters, Mq l , Mh r , M^ r , A u , Ad and thus requires 5 renormalisation conditions. For a 
physical on-shell renormalisation this requires trading these Lagrangian parameters with 5 
physical parameters. Owing to SU (2) invariance, the soft-breaking mass parameters Mq l 
of the left-chiral scalar fermions of each isospin doublet are identical. Thus, one of the 
physical squark masses, say u±, could be expressed in terms of the other masses which will 
be used as input. The mass of the U\ would then receive a finite shift at the one- loop level. 
In order to find the physical fields q^, we introduce a rotation matrix R q such as 

" 1 )=rJ R q =( <*' S9 «). (3.8) 

<?2 J \ Qr J V ~ s e q c o q J 

This transformation diagonalises the mass matrix -M~, 

Mf = R q M 2 R\ = diag^, m%) , m| > m 2 2 . (3.9) 
The physical masses are expressed in terms of the soft-susy mass terms as 



l/„o „o , . 1 



< 2 = 2 [ M Ill + M 2 RR ) ± - yf [M\ LL - M 2 - RR Y + A{M 2 LR Y ) . 

For further reference it is useful to express s 2 6» g ? the parameter that measures the amount 
of mixing, in terms of the Lagrangian mixing parameter and the physical masses 
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Note also the trivial fact that S2e q as expressed through Eq. (13.111) is regular in the limit 
— > m~ 2 since this necessarily corresponds to no mixing with M~ LR = 0. In this limit 
we can take 9 q = 0. At tree-level S2e q can be accessed directly through the q\ — > q^Z^ (or 
Z°* — > q±q 2 ) which is described by the Lagrangian 

C M2 z = igzTf-f- ([jfd~f2 + ffdh)Z°)j . (3.12) 

Provided both parameters Q u ^ have been determined along side the physical masses of 
U2,d-2, d\ one determines the tree- level U\ mass 

m li = ~T ( c e d m l + s e d m % ~ s l™\ % +m 2 u -m 2 d + c 2l3 M^) , (3.13) 

In principle we could also use all four masses as input and trade this input with one of the 
mixing parameters, leading to 

2 c l m j + s l m j ~ m l +m 2 u -m 2 d + C2pM 2 w 

m u 2 ~ m k 

Note however that the appearance of {m\ 2 — m\^) in the denominator makes this definition 
subject to large uncertainties especially for nearly degenerate masses of u\ and U2- The 
definition from a decay such as u\ — > U2Z , if open, is more direct. This is reminiscent of 
our discussion about the choice of a good definition of the parameter tan (3 in pQ . Compared 
to the case of the neutralino/chargino system, the extraction of the underlying parameters 
in terms of the physical mass parameters is rather trivial. In fact the most important 
underlying parameter to extract here is Af as this will be needed for the coupling to 
Higgses. 

3.1.2 Counterterms 

So far all fields and the parameters of the Lagrangian should be considered as bare quanti- 
ties. The bare parameters for example labeled as Vq will now be split into a renormalised 
parameter V and its counterterm 5Vo. 

It is very important to stress that the rotation matrix is defined as renormalised in our 
approach. This we have pursued consistently throughout all the sectors. Therefore from 
Eq. d£ED 
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This allows to introduce the wave function renormalisation directly on the "physical" fields 
after rotation to the mass basis. These field renormalisation constants will be chosen so 
that one gets rid of the mixing introduced by the mass shifts, at least one of these physical 
particles are on their mass shell. We therefore introduce the following counterterms 

fto = (<% + 2 SZ ij)& > ( 3 - 16 ) 
M% = M\ + 8M\. (3.17) 

The shifts on the parameters induce, 

5M? l + 5 {m\ + c 2/3 (T g 3 - Q q s 2 w )M 2 z ) 5 (m q A q ) - 5 (m q ^) 

/ _2T 3 \ \ (3.18) 

S {m q A q ) - 5 (m qf itp q j 5M? R + § (m 2 q + c 2 pQ q s 2 w M 2 ) 

After shifting the parameters and the fields, the renormalised self-energies for the squarks 
are given by 

S^(<? 2 ) = S^(g 2 ) + 6ml ~ \ 5Z W ~ m D - \ 5Z W ~ O ■ ( 3 - 19 ) 
The counterterm 5m~ is connected to the counterterm <5.M~ through the relation, 

Hij Hij u ' 

5ml. = (R^M\R^ . .. (3.20) 



3.1.3 Constraining the wave function renormalisation constants 

The residue condition at the pole for the diagonal self-energy propagator imposes 4 condi- 
tions on the diagonal wave function renormalisation constants, for q = (u, d): 

5Zl = ReT! Ml {ml), 

SZl 2 = ReZ hh (ml). (3.21) 

We impose that no mixing occurs between the two squarks gi and q\ when on-shell, con- 
straining the non-diagonal wave function renormalisation constants accordingly: 

2 

5Zf 2 = — T (ReE qiq2 (m 2 q2 ) + 5m 2 g ) , 

m q2 ~ m ii 

5Zi = 2 _ (iteE^(m|) + gm^). (3.22) 

" l qi " L Q2 

In our approach, the non-diagonal wave functions are not completely determined at this 
stage because the mixing counterterm 5m 2 - i2 appears in their definitions. It is also impor- 
tant to point out that unless <5m? is chosen judiciously these non diagonal wave functions 
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are ill-defined in the limit m? x — > m 2 ^. For further reference it is interesting to define 



m~ 2 - 

5Z ?2 A = — r (^9192 0|) + ^e s 9i9~2(^l) + 25m? i2 ) , (3.23) 



m~ 2 - m qi 



such that 



6Zf 2M = 6ZfS±6Zg. (3.24) 



Only SZf^ is now potentially singular in the limit m? — > We will come back to this 
issue when fixing a renormalisation for Srri~ 12 . 

3.1.4 Renormalisation of the mass parameters, physical masses as input 

The default scheme in SloopS takes m;, and m i2 (the lightest up-type squark) as 
input parameters considered to be the physical masses of d\, d,2 and Ui respectively. This 
fixes 3 counterterms: 



5m d n = -^^K-J 

8m l 22 = - Re ^d 2 d 2 {™\) 



6m l 22 = - Re ^u 2 u 2 (ml 2 ). (3.25) 

3.1.5 Renormalisation of the mass parameters, the issue of the mixing param- 
eter at one-loop 



To complete the renormalisation of the squark sector for each generation, as we need 5 
renormalisation conditions, we have to impose two additional conditions on what measures 
the mixing in the up squarks and down squarks and therefore fixes Sm^ i2 for tildeq = u, d. 
Once this is fixed, the remaining heaviest up squark u\ mass receives a finite correction at 
one-loop. One possibility is to define these mixing parameters through physical observables. 
One can for example choose the two decays d\ — > and d\ — > u-iW~ as inputs provided 
they are open. This is within the spirit we have followed to define a gauge-invariant tan/3 
from the decay A — > t + t~ pp. This will then define Ad and A u at one-loop respectively. 
The one-loop radiative corrections to sfermions into gauge bosons have been studied in 
previous work [TSJ, [TO]. Since the issue of mixing is quite subtle with many definitions 
based on two-point functions being rather ad-hoc, we look at the problem rather afresh. 
Moreover the discussion is the same for sleptons with mixing, we therefore generalise this 
for sfermions in general and consider that the counterterm 5m 2 ? absorbs the ultra-violet 

112 

divergence of the decay fi — > /2Z . We have just seen for example that at tree-level this 
coupling is a direct measure of the mixing. Taking a physical observable will unravel how 
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to possibly extract a gauge invariant universal definition based on the two-point functions. 
With Ai representing the tree-level amplitude, A4 = igzTj?S2o f /2, the one-loop correction 
can be written as 

M kh*> = + + 

\ e c w s w * * Z J 

(- 

+tg z Tf {l-As 2 w \Q f \) - 



/M /i/2 (mp-M /i/2 (m;- 




2 2 

\ m- — m- 

\ n h 

(25m 2 - + ReE; ?(mt)+ReEf f (m--; 

+i dZ Tf c 29f — ^ ™±J*L ) . (3.26) 

1 1 \ m- — m- 

\ h h 

The first part of the correction proportional to the tree-level contribution is due to diagonal 
wave function renormalisation and renormalisation of the gauge parameters. Just like the 
tree- level contribution this part is regular in the limit [m 2 1 — ) — > 0, see the trivial 

remark we made after Eq. (13. lip . Sy^ 2Z represents purely one- loop virtual corrections 
which do not necessarily vanish in the limit of a vanishing tree-level mixing with Of = 
much like the one- loop induced f\ — > jVy. The corrections in the third and fourth line of 
Eq. (I3.26P are due to fi <-> f 2 transitions triggered from the diagonal couplings fifiZ. 
j\/[hfo z contains pure QED corrections that can be unambiguously extracted, these con- 
tain infra-red singularities that need to be combined with the bremsstrahlung corrections. 
Subtracting these pure QED virtual corrections and the corresponding gluonic QCD cor- 
rections defines a gauge invariant, infrared safe observable that does not depend on any 
experimental cut-off on the energy of the bremsstrahlung photon or gluon. Let us define 

this observable as A4 1 . 5m 2 - defined from M. l by requiring that the one- loop 

correction, (M.^ 2Z — M.^ 2Z ), vanishes constitutes a fully gauge invariant, although 
process dependent, definition of 5m 2 - . In this definition process dependent vertex cor- 

712 

rections combine with self-energy contributions leading to a gauge independent definition. 

Eq. 13. 261 is also instructive in that it reveals how to extract a process and gauge independent 

definition of 5m 2 - . Indeed Eq. 13.261 exhibits a specific pole structure in [m 2 - — m 2 - ). The 
/i2 * v ji h 

residue of the pole must be gauge independent. Therefore considering a Laurent series of 
the amplitude in the pole (mj~ — m^)0J, a gauge and process independent definition based 
on two-point functions can be defined as 



5m p 
hi 



- lim [ ReLf ? (m 2 ? ) + ReL? ? {m\ )] = —ReY^ * (m\ , m? ) . 
2 m 2 - ^m 2 - V •ft-' 2 h ^^ 2 f 2 J f 1 ? 2 f 1 ■f 2 



m - — >m 



(3.27) 



T This is in line with the definition of the Z° mass from e + e — > through a Laurent series based 

on analyticity properties of the S'-matrix, see |17j . 
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The value at the pole ReTf- - (m 2 ~ , w?- ) is gauge-invariant and universal. All the remaining 

Jl J2 J 1 J 2 

contributions in Eq. (I3.26P are then regular in the limit (m 2 ? — m 2 - ) — > and in particular 
the contribution in the third line of Eq. (I3.26p . 

Care should be taken in defining these limits. It is useful to express m 2 ? in terms of m 2 ~ 

fl.2 }± 

2 i 2 

m - ± m t 

in order to make the dependence in the pole m 2 ~ explicit. Then S; z (m| ) = E ? ? (m 2 : , ±m| 
so that Sj i ^ 2 (m 2 ? ) is a function of these two variables. These functions should be expanded 
in m 2 j , such that 

dE' f , (m 2 - , 0) 

V hh {m) + ,±m)_) = S /i/2 (mJ + , 0) ± m)_ hh Q J+ +■■■ (3.29) 



7- 



We then have 



i?eE f - f - (m\ ) + i?eS * « (m 2 - ) i?eS * * (m| , 0) m 2 - 



m 2 - -m 2 ~ m 2 ~ 2 ^ ^+ ^~ 

,/l 72 /- 



iteE,- * (m 2 - ) - ReEf f {m 2 ~ ' 

2 2 

m~ — m 7 

h h 



ReE'jj 2 {m 2 f + ,0) + C((m 2 -J 2 ). (3.30) 



We can identify 



i?e£? ; (m 2 ? , m 2 ? ) = ReLf ? (ml , 0) . (3.31} 

/l/2 V /l' ./V /l/2 V /+' 7 V ' 



By looking at the pole structure of the amplitude it is now clear that i?eEj i ^ 2 (m 2 f ,0) is 
gauge independent. However, ReY/r^m 2 -, 0) in Eq. (I3.30p . for example, is not guaranteed 
to be gauge independent. Its gauge dependent part cancels against those contained in the 
vertex corrections. 

One should be aware not to systematically equate 

i?eS /i/2 (m 2 + ,0) = Re^ hh {p 2 = m} + ). (3.32) 

Indeed a naive replacement Rellj x j 2 {p 2 = m 2 ~ ) may still give extra contributions that are 
of order m 2 ~ . This is exactly what happens when we calculate ReH ^^(p 2 ) in a gauge 
which is not the Feynman gauge with ^w,z ^ 1- One finds that the gauge dependent 
part of the quantity ReHf^ip 2 = m 2 - ) proportional to (1 — (,w,z) ar e of order m 2 , see 
[T8l [19]. Let us mention that the choice based on ReT l j i j 2 (p 2 = m 2 ) had been advocated 
to improve the scale independence of the mixing angle [T8] . 

Note that after the renormalisation of the mixing has been set according to Eqs. fl3.27l . l3~3Tj) . 
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the last term in Eq. (13.261) contributes an ultraviolet finite part. This, on the other 
hand, is not the case of the contribution from the third line in Eq. (I3.26p . Indeed its 
ReT,' ; 7 (p 2 = m 2 ? ) might still be needed to absorb possible infinities from the vertex vir- 

tih K f+' 

tual corrections for example. 

In SloopS we work in the Feynman gauge with £w = = 1. At one-loop S.-.;, is insen- 
sitive to the non-linear gauge parameters in Eq. 12.11 We therefore obtain the same result 
for ^f.f. as in the usual linear gauge within the Feynman gauge. Therefore one can afford 
using Eq. (13.321) . Taking this into account with Eq. (13.271) and Eq. (I3.3ip . the default 
scheme in SloopS is 

5m /l2 = _jReE /l/2^ 2 = m /+) ' ( 3 - 33 ) 

To compare with results in the literature we have also implemented the prescription, 

6m\ 2 = -i(i?eS /i/2 (mJ-) + J ReS /i/2 (m|)) . (3.34) 

which is equivalent to the condition introduced in Ref. [20]. In the Feynman gauge the 
difference with the default scheme is ultraviolet safe and numerically small, see the examples 
in Section [5TT1 and [5751 

As we stressed repeatedly, in our approach we do not introduce counterterms to the rotation 
matrices since non-diagonal wave function renormalisation is necessary in any case. For the 
sfermions this reveals more easily the correct prescription to take for the renormalisation of 
the mixing parameter. In practically all other approaches counterterms to mixing matrices 
are introduced and therefore Of — > 9f + 89f. We can recover these approaches by, for 
example, looking at the example of fi — > and considering the shift to the angle, 

rather than introducing the shift 5m 2 - indirectly through the non-diagonal wave function 
renormalisation constants. From 5s2e f = 2c2e f 58f we make the identification 

5m 2 -. 

56 f = _ /l2 _ . (3.35) 
m z ~ — mi 
h h 



3.1.6 SUSY QCD corrections and the squark mixing angle 

There have been many proposals in defining this angle or alternatively the mixing param- 
eter when considering purely supersymmetric QCD corrections. The different proposals 
relied on constraining the mixing angle, Eq. I3.35[ through a combination of two-point func- 
tions in order that some specific observable be finite. This rather ad hoc approach would of 
course guarantee finiteness for that observables but does not necessarily guarantee that this 
observable or quantity is gauge invariant with this choice of counterterm. What is worse is 
that if one uses the same prescription when considering one-loop electroweak corrections 
to the same quantity even finiteness is lost. The prescription based on the residue of the 
pole would have given the correct procedure. The aim of this subsection is to understand 
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why finiteness is obtained in the case of supersymmetric QCD corrections. 
Pure QCD contributions to are from the gluino g exchange self energies and the tad- 
pole squark exchange. The results can be written in a very compact form, see for example 

m 

S L~> 2 ) = ^9 q S2e q U (ml) - A (ml)y (3.36) 

The loop functions A and B are as defined in [Hj. It is evident that E^-^p 2 ) is of order 
m? 2 — m 2 ^. It independently vanishes for s 2 e q — > 0. Note that the QCD contribution of 
the gluino does not depend on the squark masses for a general p 2 . Therefore, ^ ~^~ 2 (m 2 ^ ) — 
E~^ 2 (m~ 2 ) is finite. This explains why different schemes work fine, in the sense of leading 
to finite results, for SUSY QCD corrections to processes involving squarks. One of the 
most complicated is based on tuning combinations of E^ such that a finite results for 
e + e~ — > qiq 2 obtains as far as QCD corrections are concerned [22]. With the coupling of 
the Z to squarks defined as Cy for Zq^qj, the following combination is used to define the 
counterterm, 



C22 — Cn 



This can be rewritten as 

J qiq 2 ( m gi) — c n j 



c 22 Re^ 1 a 2 (ml ) - cnReE&s* (ml ) ReE*, & (m? ) + EeEg, „- 2 (m\ 



(3.37) 



(3.38) 



C22 - c u 2 

c 22 + c n geggigjwjj ~ ReEq^m 2 ^) 
C22 - cn 2 

The much simpler scheme based on the use of i?eEq 1 g 2 (m?) [23] is in fact a very special 
case of the scheme in Eq. (I3.38p . we can see that it is obtained as Cn — > in Eq. (13.381) . 
For the electroweak case the extra terms proportional to ReH^^mi-J — -ReEg 1 q 2 (m? 2 ) in 
Eq. (I3.38P are not finite apart from the gauge invariance issue. However as we have seen the 
ultraviolet divergent part can be cancelled in (i?eEq 1 g 2 (m| 1 )-|-i?eEq 1 q 2 (m~ 2 ))/2 as suggested 
in [21]. However this suggestion was not based on a very strong theoretical or physical 
argument apart from it being more symmetric or democratic in the two squarks. 



3.1.7 Deriving the counterterms 

We are now in a position to derive all the needed counterterms. First of all with both 
prescriptions for Sm~ l2 either based on Eq. (13.331) or the naive Eq. (I3.34[) . the non-diagonal 

wave function renormalisation constants 5Z% and 5Zf, are now regular in the limit m? — > 
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m 



| , where any potential ultraviolet divergence is contained in 5Zf 2 ■ In fact in the scheme 



q,S 
12 • 



of Eq. (13.341) only this part remains and therefore 5Zf 2 = 5Z 21 = ^ 
The remaining counterterm Sm^ is completely constrained, 

6m L = -IT fe m L + s ejrn 2 d - s 2 e d 5m\ 2 - s 2 g Jml 22 + s 29u 5ml i2 

+ 8ml - + M 2 W (cj-^f - s 2 Jj^) ) , (3.39) 

For Cg 1 this scheme is not appropriate as it will induce large radiative corrections. 
One should prefer the use of m\ x as input parameter in lieu of m| 2 . With m\ 2 as input 
parameter the physical mass of u\ will then receive a finite correction at one-loop, 

= ml + 5ml n + iteE^ (t< ) . (3.40) 

Alternatively we can use m„ 1 as input like we have done with the other squark masses. 
This will allow to define tp from the sfermion sector through 



— - = o ; r0 (ci 5m 2 , + si 5m 2 r — s 2 a,5m~ —5m 
in s M V d dl1 d d22 dl2 



5tt 

t/3 £>2/3 IVI W 

- c e u 5m ln + s lj m l 22 ~ s 2e u 5ml 12 + 5m 2 u + c 2(3 5M^j . (3.41) 

Using Eq. ( 13.201) . we find the relations between the counterterms Srn^ and the countert- 
erms 5Mq l , 5Mu a , 5M^ r , 5A u and 5 Ad of the underlying parameters at the Lagrangian 
level 



SM. 



(ci 5m 2 . + si 5m~ — s 2 a.5m\ — 5m 2 , 

Ql 

6MilR = 2M~ ( s «« 5m ^i + $v Sm *22 + s 2fl u 5m? 12 - 5m 2 t 

2 2 A/ r2 ( n ( 5M Z , Ss w\ 2 8tp 



5Mi = ( si 5m 2 j + ci 5m 2 r + s 2 Q,6m~ — 2m d 5m 

R 2M- V 11 22 12 



5{m u A u ) = —(5m iLii -5m iL22 )+c 2du 5m ili2 + —\5{i + ii y,-f- , 

5(m d A d ) = "-f (5m\ u - 5m\^ + c 2dd 5m\ 2 + ™ d tp (sp + ^ + ^ (3.42) 
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3.2 Slepton sector 



After having shown in detail how the squark sector is renormalised in the case of mixing, it 
is straightforward to treat the slepton sector. Again for the sleptons, the case with mixing 
is for all practical purposes only relevant for the f . In the code we have implemented 
mixing for all generations, in the first and second generation this is used only in to conduct 
high precision checks on the results, for applications the unmixed case is used. Here we 
will show only the case with mixing, the unmixed case is then trivial. 
Compared to the squark sector, as seen from Eq. (13. 2p . one has, for each generation, only 3 
parameters : Mi , Mg R A e and one field is missing, u R . e L and will mix leading to the 
physical fields e\ and e 2 . In the unmixed case we associate e.\ with e^. The mixing matrix is 
described in exactly the same way as in Eq. (13.41) with q — »• e and the different components 
given by Eqs. (13.51) - (13. 7p with Q — > L with the corresponding quantum charges. Shifting 
the fields and parameters we can write the self-energies (in the case of £ diagonal and 
non-diagonal) as 

S^(g 2 ) = S^(g 2 ) + ^-^4(g 2 -m|)-^4(g 2 -m|), 
%{q 2 ) = ^(q^ + Smj-SZ^q 2 -mj). (3.43) 

We take the physical selectron masses as input parameters through the usual on-shell 
condition. We require the residue of the propagators of ej and v to be equal to unity and 
no mixing between e\ and e 2 when these are on-shell. These conditions imply 

5mj = — i?eSe e (^ ) , 

en e l e l \ ei/ i 

6Z* = iteE^mj,), 
2 

5Z 12 = ml (^(mD+fay, 
/ / L~ 1 1 L ~ 

e2 ei 

5Zl, = 2 _ (ReZ SlS2 (ml) + 5ml 12 ) . (3.44) 

" L e l " L e 2 

The remaining parameter 5m| l2 describing mixing is fixed analogously as in the squark 
sector. The default scheme in SloopS is 

6ml i2 = -ReE &lSa ((ml + m e 2 J/2) . (3.45) 

As in the squark sector, a better definition would be to relate this counterterm to a physical 
observable like the slepton decay e\ — > ^Z for example, see (I3.26p . The naive scheme 

Smj i2 = -1 (jfeE^O + ReE ilg2 (mD) , (3.46) 

is also implemented. Another possible scheme uses the mass of the sneutrino as an input 
parameter such that 

dm 2 = -ReEo(ml), (3.47) 
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and the counterterm 5m\ is given by, 



= ^ ($Mn + - ^ " Sml + M 2 W (cj-^f - ) ) . (3.48) 



However this definition is to be avoided since the mixing in the slepton sector is usually 
very small, S2e e ~ 0, even for r's which would lead to large corrections. 

The extraction of the counterterms of the parameters at the Lagrangian follows 

SMf = [ c 2 , 5m 2 - + si 5m 2 - — s 2 e 5m 2 - — 5m 2 

L L \ 



w 1 



SMgR = 2M~ I Sf) e Sm ^ + C 8e 6m e22 + ^JlTl^ ~ 5l7l e 



'5M 2 Z 5s 2 w \ 2 Sty 



S(m e A e ) = f£ (5ml, ~ Sm L) + WML + "Udp + ^ + ^ V (3.49) 

If the sneutrino mass is not used as input then it is predicted with a finite correction from 
its value at tree-level. 



m- 



r°° P 2 = < CC 2 + 2 ) - ReZ SlSl «)) + M 2 W (cj-^f - s 2 J-^j 

+ 4 e (fteZg^imlJ - ifeEg 2 e 2 (m? 2 )) - s 2 g e 5m? i2 - 5m\ . (3.50) 

In the limit of massless fermions, the term in the second line vanishes and we identify, as 
said earlier, §1 with e^. This is a very good limit for the selectron and smuon sector but 
we have to consistently take the electron and muon Yukawa couplings to zero. 



4 The chargino/neutralino sector and its renormali- 
sation 

4.1 Fields and parameters 

The charginos and neutralinos are mixtures of the spin-1/2 fermions which are part, on 
the one hand, of the two Higgses chiral multiplets, Hi^ which constitute the higgsinos, 
and on the other hand, the electroweak gauginos within the gauge supermultiplet for the 
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U(l) and SU{2) gauge groups of the Standard Model. In terms of the two-component 
(left-handed) Weyl spinors the two higgsino doublets in accordance with our definition in 
the Higgs sector P are H\ = (H®, H{) and H 2 = (H 2 , We denote the U(l) gaugino 
(bino) as B° and the SU(2) one as W\i = 0,1,2 with W ± = \{W X T iW 2 ). Due to 
electroweak symmetry breaking the electrically charged components will mix and lead to 
the charginos that will be collected as Dirac spinors X121 while the electrically neutral ones 
will mix leading to the neutralinos that will be described as Majorana fermions, Xi 2 3 4- m 
this sector soft masses enter only through the soft masses of the gauginos 

cf oft = -i (m.B'B + M 2 W l w\ , (4.1) 

which is the only source of mass for the gauginos before electroweak symmetry breaking. 
The higgsinos get a mass from the supersymmetry preserving \x term in the superpotential 

C* =lie ij H{Hi + h.c. (4.2) 

Supersymmetric gauge matter interactions lead to mass mixing terms between these states 
after symmetry breaking through 

= [hKqW^ - g'B°)Hi + Hl(gWV + g'B°)H 2 + h.c.) , (4.3) 

with r % the Pauli matrices. At this point let us give our convention on the sign of the pa- 
rameters /j,, Mi, M 2 . We can always take M 2 > since any other phase can be transformed 
away by a field redefinition in Eq. (14.11) . however because of the mixing term in Eq. (14. 3[) 
we loose the freedom to redefine the phases of the Higgsino and bino fields and hence the 
sign of \i and M\. 

The kinetic term in terms of the current fields writes as 

£kin = iW*a»(d^W l ) + iEPa^d^B ) + tH 1 a fi d, l H 1 + iH 2 o»d»E 2 . (4.4) 
Collecting all terms in the chargino mass matrix and defining 

W~\ , c ( W+ 



r R = [ i- ), r L = [ ~ + ), (4.5) 



leads to 

C c = i[^d^ R + - WrX^l + TlX^r] • (4.6) 

1 stands for the transpose operation and the mass mixing matrix is given by 

X - ( * . ) • (4.7) 



V2M w Cf3 /j 
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The system can be diagonalised by two unitary matrices U and V that define the physical 
(Weyl) fields as 



(4.8) 



In the case of CP conservation that we will cover here, we can take both U and V real. 
We write the diagonalised mass matrix X 



m-± 

x = UXV t = x l = vx l u = I *! 

m-± 

A2 



(4.9) 



m,-± are the (positive) eigenvalues of the hermitian matrix XX* with m-± < m-±. In our 

Xl,2 v ' Xl A2 

implementation in order to have positive eigenvalues we take 

detU = +1 and det^ = sign(detX) = e M with detX = M 2 /i - M^s 2 p. (4.10) 

The physical masses are also defined from the invariant basis independent quantities that 
are the trace and the determinant of the square matrice XX 1 , which give 



Xi x 2 

2 2 

m~±m~± 
Xi x 2 



(detXf 



(4.11) 



and 



m~± ~± 
Xi >x 2 



Ml + /i 2 + 2M^ 



T [(M 2 -/i 2 ) 2 + 4M^c 2 /3 + 4M^(M 2 + /i 2 + 2/iM 2S2/3 )]2 ) . (4.12) 



The corresponding chargino Dirac spinor \i (i = 1, 2) is constructed as 



X. i 



XI; 



xt= ( ) =XT* i = l,2 



(4.13) 



Similarly the Lagrangian for neutralinos writes 



(4.14) 



where 



(4.15) 
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The mass matrix Y 



Y 



/Mi 

M 2 

-M z s w C(3 M z c w cp 

\ M z s w S(3 -M z c w sp 



-MzSwC/3 M z s w s p \ 

M z c w cp -M z c w Sfs 

-n 

-IM / 



(4.16) 



can be diagonalized by an unitary complex matrix with the physical states being 
We will refer to the diagonal matrix as 



Y = N*YN^ = diag( 



Xl ' X2 7 A3' A4 y 



< m v o < m v o < m v o < m^o 

Xl X2 X3 X4 



(4.17) 



(4.18) 



Note that N can be written as JN where N is real and J = diag(ji, J3, J4). diago- 
nalises F but leads to masses that are not necessarily positive. A positive mass obtained 
with corresponds to ji — 1, a negative mass corresponds to ji = i. The corresponding 
neutralino (4-component) Majorana spinor Xi (i — 1)2,3,4) is given by 



v" 
A I 

v! 1 
A / 



(4.19) 



4.2 Renormalisation: Counterterms and Self-energies 

We could have treated the chargino and neutralino system that we have just exposed within 
a common notation, deriving in a compact form the neutralino sector on the basis of its 
Majorana nature. This could have been done through a mass matrix M that stands for 
either X (charginos) or Y of the neutralinos and the two fields ipR ^ that represent either 
ijjji L or the single Majorana field ip n . To make our renormalisation procedure of this sector 
as transparent as possible we will take this common approach to show that the approach in 
renormalising the chargino and neutralino sector is exactly the same and that it corresponds 
to the approach that we have taken in the Higgs sector and the sfermion sector as concerns 
the issue of mixing. In particular we stress that we do not renormalise the rotation matrices 
that express the mass eigenstates from the current eigenstates. Summarising what we have 
just seen in the sfermion sector and splitting as usual the bare Lagrangian (denoted by ) 
into a renormalised Lagrangian and counterterms, the kinetic term and the mass term of 
a fermion field ip with an arbitrary number of components can be written as 

4 = iWk ^d^ R + tf L a^ L 0] - W R M ij R + j? L M$ R ] , (4.20) 

where i/jr/l represents the the fermion field and Mq the non-diagonal mass matrix at bare 
level. At tree-level this mass matrix is diagonalised by rotating the fields with two unitary 



20 



matrices D R and Dl which define the current fields so that at bare level we write these 
fields as 

XR0 = D R ip R0 , Xlo = D l ^ lo . (4.21) 

The corresponding diagonal mass matrix M is then given by, 

M = D* R MD[ = M+ = D L M l D R = diag(m^,m^ 2 , ...) , (4.22) 

and gives the physical masses m^. The ensuing Dirac/Majorana spinors x%o are con- 
structed with these Weyl spinors 

X,o = if') • (4-23) 

After the diagonalisation is performed, the counterterms for the different parameters in- 
volved in the mass matrix are introduced, 

M = M + 5M, (4.24) 

and also the wave function renormalisation constants 5Z-' for each chiral "physical" field 

Xr/l, 

XR,Li\ = (Sij + \ 5z 5 ,L ^XR,Lj- (4.25) 

These transformations for the chiral fields are equivalent to the following transformation 
for the four-component spinor x%i 

X l0 = X* + \ [SZ^Pl + 6Z%*P R ] xj . (4.26) 

We stress again that in our renormalisation scheme, we do not use the extra shifts on the 
diagonalisation matrices D L R , D L R — * D L R + 5D L R , in other words 5D L R = as done 
in Ref. [25]. This is in the same spirit as within the Higgs sector and the sfermion sector. 
So, we consider that the diagonalisation matrices D^ R at tree-level and at the one-loop 
level are the same, Dl R are renormalised. With the renormalisation counterterms (14.241) . 
fl435|) and 

5M = D* R 5M D\ , (4.27) 
the renormalised self energies S^iXi can ^ e cas ^ 

>W?) = ^xA<l) - PlSM.j - P^M*, (4.28) 

- m^[5Z^P L + 8Z«*P R ] + \[8Zf*P R + 5Z*P L ](4 - m % ) . 
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Eq. 14.281 shows clearly that the wave function renormalisation constants are not involved 
in the renormalisation of the Lagrangian parameters contained in the mass matrices M 
which in our case involve Mi, M2, [i. 

It is useful to decompose the self-energy into the independent Lorentz structures through 
the projectors Pl,r = ^f 52 '; 

e** (?) = P£^u ^ + p * s & ^ + i p ^iu ^ + ^ ■ ( 4 - 29 ) 

Hermiticity imposes the following constraints on the elements of the Lorentz decomposition 

These are also satisfied by the corresponding covariants of the renormalised self-energies 
in Eq. 14.281 For a Majorana fermion (like a neutralino in the following), the additional 
Majorana symmetry imposes 

*fttf> = E&tf). = E&tf). q&tf) = = s^tf) ■("•31) 

Some of these properties are used in our code as an extra test. 

To fix the wave function renormalisation constants 5Z^' L , we require that 

• the propagators of all the charginos and neutralinos are properly normalised with 
residue of 1 at the pole mass. This pole mass may get one-loop correction. For our 
treatment at one-loop it is sufficient to impose the residue condition by taking the 
tree-level mass. Taking the one-loop mass is a higher order effect, see Section 4.7 of 
[T] of our treatment in the Higgs sector. This condition implies 



,ml t q 2 -m x 



, lim 2 ^ — df^XiXi (?)«*(?) = u xM and 



lim u x MR^mM^§- = (4-32) 

H Xi ^ Xi 

• No mixing between the physical fields when these are on mass-shell 

Ret^ 3 {q)u Xj (q) = for q 2 = mj. , (z ^ j) . (4.33) 

With these conditions wejio not have to consider any loop correction on the external legs. 
Note that as usual [H] Re signifies that the imaginary dispersive part of the loop function 
is discarded so as to maintain hermiticity at one- loop. Eq. 14.331 gives the diagonal element 
of the wave function renormalisation constants 



SZff = -ReEfC (ml )-m%( ReZ± v ~ (m| ) + ReHf- (m% )) - 2m y ReL 

« XiXi y Xi' Xi V XiXi y Xi ' XiXi y Xi' j Xi 



? s - (m 2 - 

XiXiV n Xi 



(4.34) 
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where we have used the fact that in the case of CP conservation im% ) = Ti RS ^ (m~ 

XiXi v Xi ' XiXi v Xi ■ 

The prime on a function such as E|^.(m|.) stands for the derivative dT, R ^.(q 2 ) /dq 2 \ q 2 =m 2 
The non diagonal elements (i ^ j) of 5Z L,R are derived from the constraints of Eq. 14.331 

SZ^ = — = =- [rn* ( R&Tk% (m% ) + my.ReL R %.(m% ) + m^rriy . ReE RX ~ (ml 

*J m 2 _ ^2 ^ Xi XiXj\ Xj' Xj XiXi v Xj/ Xi X] XiXj\ X 



Xi Xj 

+ m\. ReZ^ ( m % )~ m xr 6 Mij - m Xj 5 M* { . 



SZ t = m 2 2 _ m 2 (™x^ S &K) + m Xi ReE R ^(m%) + w^ReE™,^) 

Xi "Xj 

+ ^^^y^)-^.^-^^)- ( 4 - 35 ) 



Specialising to the case of the charginos we will have to take Dr = U and Dl = V, see 
Eq. (14. 7p and M — X, see Eq. ( 14. 8 p where both U and V are real matrices as is the mass 
matrix X in our case with CP conservation. In this case 5Z L,R can be taken real. For 
the neutralinos Dl = Dr = N, see Eq. ( I4.17P and M = Y is a symmetric real matrix, see 
Eq. KIB . In this expected, we have 6Z L = 6Z R = 6Z° which is a result of the 

symmetry of Y and the Majorana constraints of Eq. 14.311 In fact Eq. 14.351 can be recast 
into 

= (m. x o(Re&(m 2 ) + Re^t(m 2 

Xi Xj v 

+ (^£go(m|o) + ^g(m|o)) - (SYij + 8Y* 
+ ^ ( -m x0 (ReE^Am 2 ) - R^ L lA{m\ )) 

Xi Xj \ 

+ (SeE^^o) - feSgK;)) - - ^)). (4.36) 

Note that while Z\ is real, Zfj, (i ^ j) is either purely real (given by the first two lines in 
Eq. 14.361) or purely imaginary (given by the last two lines in Eq. 14.361) when using N = JN 
in order to have positive masses. Using N = N we can have Zfj real but we have to keep 
track of the sign of the masses in Eq. 14.361 For example, with N = N (real), when both 
m Xi . > are obtained, we get 

&Z% = ( ^5^(4) + ^^^(mL) - ^y) . (4.37) 

m o — 777,- V Xi Xj Xj x 3 Xi Xj^ Xj' J 

Xi Xj N / 

It is important to note a common feature of our approach that we already encountered in 
the case of the mixing in the Higgs sector and the sfermion sector. The non diagonal wave 
function renormalisation constants in Eqs. I4.35[ I4.36[ 14.371 are fully determined only once 
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the mass counterterms SM are fixed. In our case this requires fixing 8 Mi, 8M2 and 5/x to 
which we turn in the next section. 

For completeness let us give the corresponding counterterm matrices SX and SY. We have 



SX 



with 



8M 2 8X U 
5X 2 i 5fi 



5X U 
5X 2 i 
8Y 13 -- 

5Y 14 - 
8Y 23 ~- 
5Y 24 -- 



8Y 



5Mi 





SY 13 


5Y U 





6M 2 


SY 23 


8Y 24 


8Y t3 


8Y 23 





—Sfi 


SY 14 


SY 24 


—6/1 






±V2M wC/3 
-s w M z c p 
-swMzSp 

-c w m z c p r 

-c w M z s p 





1 Sslf 


Q 2 5t 






1 <5s»- 




1 Sc%^ 


„2 % 


2 ^ 




1 $Cw 
2 ^ 


+4t 



(4.38) 



(4.39) 



4.3 Fixing 5M h 5M 2 ,5fi 

5Mi,8M2,5fi can be fixed through the diagonal self-energies of the chargino-neutralino 
system which we have not fully exploited yet and which constrain the physical masses 
of the charginos and neutralinos. The most straightforward and simple choice is based 
on the fact that the chargino system is a 2 x 2 system which is easier to handle that 
the 4x4 system of the neutralinos. In SloopS the default scheme is to choose the two 
chargino masses m~± and m~± as inputs to define the two parameters M 2 and u and one 

Xi X2 

neutralino mass to define M\. The lightest neutralino mass m^o is used by default to fix 
Mi. The three other neutralino masses m^o are derived and receive one-loop quantum 
corrections. At one-loop these three input parameters translate into the usual definition 
of the pole masses in the on-shell scheme through the renormalised self-energies of the 
charginos and the lightest neutralino 

Re ^xixM) u xM) = 0for 9 2 = m xi' for X* ~* Xi,xf,X°v (4-40) 
This translates into 

11 Xi xfxi xf 2 Xl xfxf v xf y xfxf v xf JJ 

6X22 = 8^ = M^(m^) + ir^(fe^(m^) + fe^(m^)), 
8Y n = 5 m ,o = ^o(my+m^eEgoK„). (4.41) 
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These three counterterms can be inverted to derive the counterterms parameters 5 Mi, 5M2, 
5n through 5Y = N*5YN^ and 5X = U*5XV\ see Eq. ( |4~2T1) . In fact SM 2 , 6/1 can be 
derived more directly without going through the mixing matrices from Eq. (14.121) . We get 

1 / Sm~+ 5m~+ 

SM 2 = - (M 2 m 2 + - udetX) — ^ + (M 2 m 2 + - /idetX) — ^ 

M, — a 2 \ *i m.-+ x 2 m + 

\ Xi X.2 



M w tp J 



5m~+ n 5m-+ 

5 V = 2 \,2 I (V m i+ ~ M 2 detX)^± + (fim 2 + - M 2 detX)- 



u? — M, 2 I *i m~+ *a - m 



xi x 2 
M^(/i + M 2 s 2f3 ) 6 -^f - M 2 M 2 w s 2 pcJ-^\ , (4.42) 
1 



= W2 (Sm x o-N* 2 2 5M 2 + 2N* 13 N* u 5v 



ii 



N- 

-2N* n N* 3 5Yi 3 - 2N* 2 N* 3 5Y 23 - 2N* n N* u 6Y u - 2N* 2 N* U 5Y 24 ) . (4.43) 

The physical masses of the other three neutralinos (i = 2,3,4) receive a correction at 
one-loop given by 

m^o ys = rriyo + 5Yu — ReY,%o- (rr& ) — m^qReY^~ {rn%a) with 

Xi A a Xi Xi Xi XiXi Xi 

6Y tt = N* l 2 5Mi + N* 2 2 8M 2 -2N* 3 N* A 8fi 

+2N; 1 N* 3 5Y 13 + 2N* 1 N* i 5Y u + 2N* 2 N* 3 5Y 23 + 2N* 3 N* A 8Y 2A . (4.44) 

Checking the cancellation of the ultraviolet divergences in Eq. (I4.44p is an important non 
trivial test on the validity and correctness of the procedure and its implementation. 
Other schemes in the neutralino/chargino sector can be implemented in SloopS as will be 
shown in a forthcoming publication. A deviation from the commonly used scheme adopted 
here was taken in Ref. [26] where the input parameters are the masses of x\i X2 an d xt- 
There are a few important remarks to make about Eq. (14.421) and Eq. (14. 43 p . The choice of 
m^o as an input parameter is appropriate only if the lightest neutralino is mostly bino or 
if the bino-like neutralino is not too heavy compared to the other neutralinos. Otherwise 
the extraction of Mi would be subject to uncertainties. This shows in Eq. (14.431) since 
iVii would be too small which would in turn induce large radiative corrections. Another 
difficulty arises with the special configuration M 2 ~ ±/i. Eq. (14.421) shows that an apparent 
singularity might be present. We had already pointed out in [6] that this configuration 
can induce a large tg-scheme dependence in the counterterms 5M 12 and S/i and therefore 
to the annihilation of the LSP into W's for a mixed LSP, see also H7\ . Let us look at this 
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configuration again. We can rewrite Eq. (I4.42p as 



$M 2 = j^—^ [e^5E x + (M 2 - e^)5F x j = -^—^ [\^E X + (M 2 - \p\)6F x ) , 
Sfi = [e,M 2 8E x + (ji - e,M 2 )5F x ^} = (m 2 8E x + (|/i| - M 2 )5F X ^ , 



5F X 



i (5ml t + 5ml t ) - 5M^ . (4.45) 



It is important to note that the contributions proportional to 5F X are regular in the limit 
M 2 — > \fj,\, moreover 5F X does not introduce any tp dependence. Only terms in 5E X may 
cause trouble. The problem is confined to the finite part (in the ultraviolet sense) of SE X . 
Indeed, we have checked explicitly that in the limit M 2 — > 5E X is finite. This is a strong 
check on the validity of the code. Therefore any non regular term comes from the finite 
part (in the ultraviolet sense) of SE X and calls for a good choice of the renormalisation 
scheme in order not to induce too large corrections or ill-defined constants. 



4.4 Input parameters and parameter reconstruction 

In practise, in the on-shell scheme that is generally used for the chargino/neutralino sector 
and that we adopt here we need to reconstruct from experiments the value of /x, M 2 and 
Mi from three physical masses. If we invert the mass relations of the chargino sector, we 
would in general get four solutions (M 2 ,/i) for one set of chargino masses (m-+,m-+) 

1 1/2 



m\ + m\ - 2M^ 



ml+ + m|+ - 2M^) 2 - A{M^s 2 p + e M m-+r? 



2 Xj X 2 _ 

2 2 l K ""xZ 1 ""xl —wj -\—w-w ' ~n~~xi"~X2 

M 2 = [ml + +m\ + -2M 2 w -^] l l\ (4.46) 

where e MiX can take the value ±1 and summarize the ambiguities in the reconstruction [28] . 
e M represents the sign of [i so that [i = e^\/~jj?. e x represents the M 2 <-> /j symmetry in the 
reconstruction so that Sgne x = Sgn(/i 2 — Mf). In the numerical computations of the one- 
loop correction to the neutralino masses in Section l5\2l we have taken the set corresponding 
to e x = e M = 1. Once M 2 and \i are known, the remaining parameter M\ can be extracted 
from the knowledge of one of the masses of the neutralinos. For example, in the case where 
the neutralino is mostly bino-like and corresponds to the lightest neutralino with mass m x o 
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as what occurs with the models with gaugino mass unification at the GUT scale, we have@ 

mi, - M 2 m 3 - (/i 2 + Mf Wo - (s 2/3 Mf /i - (/i 2 + s 2 w M 2 z )M 2 Wo + s 2p s 2 w M 2 ^M 2 
1 " m| - M 2 m 2 - (/i 2 + c 2 w M 2 )m^ - s w c 2 w M 2 ^ + /i 2 M 2 

Al A-l ' t l 

(4.47) 

Having Mi, M 2 , fj,, a consistency check can be made to make sure that M\ is indeed given 
through Eq. (I4.47P with m x o as input and not some other neutralino. This shows somehow 
the ambiguity, already encountered in extracting M 2 , \i from the 2 chargino masses, in 
reconstructing the Lagrangian parameters from the knowledge of three masses only. This 
said, considering that, with the present limits on the chargino masses, the effect of mass 
splitting is small like, as we will see, the effect of the radiative corrections on the neutralino 
masses, discovery of both charginos almost certainly guarantees the discovery of the two 
Higgsino and the wino-like neutralinos with masses of the same order as the corresponding 
charginos, therefore allowing to select the correct (M 2 , fx) from the chargino reconstruction. 
If the bino like is not too heavy it will then be easy to single out and hence measure M\. 
Another exploration about the correct extraction of M 2 , //, Mi can also be done through 
the measurements of some couplings of the charginos (see for example [29] for a tree-level 
analysis) and the neutralinos (see for example [3D])- We will see below how one can extract 
these parameters in decays involving the neutralinos combined with the measurements of 
the chargino masses. Although the situation here is quite different from the mixing in the 
sfermions, exploiting decays as inputs, to fix the underlying parameters less unambiguously 
when mixing takes place is promising. We will get back to this issue in a forthcoming 
publication. Meanwhile let us give an example about the reconstruction. As an example 
the measured masses that we take as input are the two chargino masses with m~+ = 232 
GeV and m.-+ = 426 GeV and the lightest neutralino mass m^o = 98 GeV. From the 
chargino masses we obtain four solutions for (M 2 ,/x) according to Eq. (14.461) . For each one 
of these solutions we first reconstruct the corresponding Mi by imposing the mass of the 
Xi taken as input. In the example we have taken, the four solutions for (M 2 , n, Mi) are 
(all given in GeV) 

(250.39, 399.78, 100.38) for e x = 1, e M = 1 , 
(240.39, -405.86, 98.22) for e x = 1, e M = -1 , 
(399.78, 250.39, 103.68) for e x = -1, e M = 1 , 

(405.86, -240.39, 100.05) for e x = -1, e M = -1 . (4.48) 

Each solution will lead to different predictions on the observables in the chargino and neu- 
tralino sector as well as the sfermion/Higgs sector. Comparing the theoretical predictions 
to the measurements of a minimal set of these observables lifts the four-fold ambiguity. 
Theoretically with each set of solutions in Eq. (I4.48P we can give one-loop predictions 

•^Equation 14.471 was derived in [25j however there is a typo, in the last term in the numerator of 
Eq. ((4747)) is missing in [25] . 
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given some other model parameters that indirectly enter in the one-loop calculation. For 
simplicity and to avoid having to deal with QED corrections we consider the prediction 
on the 3 other neutralino masses and the decays X2 ~ ^ Xi(7>^°)- The former is a pure 
one-loop effect. We take the pseudo-scalar mass M^o = 300 GeV, a common soft-susy 
sfermion mass Mj = 500 GeV, a common Af = 0, the SU(3) gaugino mass is set at 
M 3 = fOOO GeV and tp = 10. For t@ the results we present below are within the MH 
scheme. <5r(x2 — ► Xi^°) is the one-loop correction to the rate X2 ~~ * Xi^°- 

The results in Table 14^1 show that disentangling between the possible solutions is in princi- 
ple possible even if not all neutralino masses are measured. For example the rate X2 — > Xi^ 
is a clear cut indicator for the sign of e x , since this rate is an order of magnitude larger 
if the higgsino-like neutralino is lighter than the wino-like neutralino. If a precision mea- 
surement below the 10% level can be achieved on this observable it can, by itself, also 
disentangle between all four solutions. Considering the smallness of the rate \X —> Xi^° 
this observable is perhaps of academic interest. Note however that it can in principle be 
used to lift the degeneracy between all four solutions. Combining measurements like this 
with measurements of some of the other neutralino masses or measuring all the neutralino 
masses is certainly a good way to lift the ambiguity. 



5 Applications and examples at one-loop 

Our code has been checked extensively. We have written a script that automatically cal- 
culates cross sections for all 2 — > 2 process in the MSSM at one-loop. We check ultraviolet 
finiteness as well as the independence in each of the non-linear gauge parameters. Results 
of these extensive checks can be found in [HT] . 

Moreover we have compared the results of the code and the renormalisation procedure with 
quite a few observables that have appeared in the literature. Apart from these comparisons 
which we will report here the flexibility of the code allows us to study the scheme depen- 
dence of the result. We show here a few examples, taken from studies by different groups, 
of comparisons ranging from mass corrections, two-body decays as well as 2 — > 2 processes 
paying a particular attention to the important tp scheme dependence. For the latter we 
consider the schemes introduced in [1] and summarised in Section [2j The examples we will 
review here cover the sectors we studied in this paper, leaving aside the Higgs sector that 
we studied at length in [Tj. 

Before embarking on showing our results for some observables at one-loop, let us briefly 
describe how we treat infrared divergences. The one-loop corrections can still contain 
infrared divergences due to photon virtual exchanges. These are regulated by a small 
photon mass. The photon mass regulator contribution contained in the virtual correction 
should cancel exactly against the one present in the photon final state radiation. The 
photonic contribution is in fact split into a soft part, where the photon energy is less 
than some small cut-off k CJ M s °f t (E 1 < k c ) and a hard part with M.^ ard {E 1 > k c ). The 
former requires a photon mass regulator. We use the usual universal factorised form with a 
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e x = !> 6 M 



tree-level = 232.34,777^ = 232.19 

X2 X2 

m tree-level = 40 5.26 , TB 1 ^ = 407.41 

*3 X3 

tree-level = 425.69,777^= 425.77. 

X4 ' X4 



r(x°^X?7) = 0.308 xlO" 8 , 
r(x° -> x?Z°) tree " level = 0.223 x 10~ 2 , ST(x° 2 -> = 0.533 x 10" 4 



e x = !> 



m ge-ievei = 2 3 1 .83 , m p | ys = 231 .74 

m t_ree-level = 4 1 4 . 2 , 777 P | yS = 414.19 
^tree-level 



422.79, m P o ys = 423.46 

X 4 ' X4 

r(X2-^X?7) = 0.182 xHT 7 
r(x° -> x?^°) tree " level = 0.202 x 10" 2 , 5r(x° -> = 0.780 x 10" 4 



e x = -1, e M = -1 



^tree-level 



: 236.17, m P o ys = 236.17, 

X2 X2 

tree-level = 256.54,777^ = 254.71, 

A3 X3 



m tree-level = 425.00,777^= 425.81, 

A.4 X4 

T(X2^X?7) = 0.142 xlO" 7 , 

r(x° -> x?^°) trce " level = 0.197 x io" 1 , 5r(x° -> = 0.271 x 10- 



mi 



^tree-level 
^tree-level 

m-o 

^tree-level 
ill -0 

x 4 



nxl - xh) 

~0 ryO \ tree-level 



231.76,777 P o ys = 232.59 

' X2 

249.64,777 P o ys 



= 249.53 , 
= 425.65, 
0.368 x 10" 10 , 

0.277 x Hr\5r(x°->x?Z°) 



X3 

425.80, m p ^ ys 
X4 



0.157 x 10 



-2 



Table 1: Disentangling between the four solutions for M 2 , /i and Mi from the input with 
777-+ = 232 GeV, m~+ = 426 GeV and m^o = 98 GeV. All masses and decay widths are in 

Xi ' Xi Xi " 

GeV units. ST is the one-loop correction in the MH -scheme. 
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simple rescaling for the case of the gluon correction in all processes we have studied where 
the non-abelian coupling of the gluon is not at play. The test on the infrared fmiteness 
is performed by including both the loop and the soft bremsstrahlung contributions and 
checking that there is no dependence on the fictitious photon mass A 7 or gluon mass A g . 
For the bremsstrahlung part we use VEGAS adaptive Monte Carlo integration package 
provided in the FFL bundle and verify the result of the cross section against CompHep [32] . 
We choose k c small enough and check the stability and independence of the result with 
respect to k c . 

5.1 Corrections to the sbottom and stau masses 

We compare our results with those of Ref. [33] where an approach similar to ours in this 
sector is taken. For tp, the authors of [33] take a DCPR scheme and compare with DR. 
The mixing parameter in [33] is however defined through the naive scheme of Eq. (13.341) . 
In order to conduct this comparison we first need to implement the same set of input 
parameters as in [33]. We therefore slightly change our scheme to predict the heaviest 
sbottom mass m~ b at one-loop instead of the heaviest stop mass which is therefore 
taken as input. Our code being quite flexible this change can be made very easily. The set 
of parameters corresponds to the (tree-level) choice /i = 100 GeV, M 1 = 95 GeV, M 2 = 200 
GeV, M 3 = 719 GeV, M A o = 150 GeV and Mf R = M f ~ L = A f = 300 GeV. This assumes 
implicitly that these Lagrangian parameters have been reconstructed from the physical 
inputs. Let us discuss our results first, taking the same scheme for the sfermion mixing 
parameter as in Ref. [33] before commenting on the impact of taking the SloopS default 
scheme for this parameter. As Fig. [1] shows, the corrections are almost insensitive to the 
ig-scheme in the case of the correction to the sbottom mass, which is very welcome. Indeed, 
the A TT , DR and DCPR are within 0.03% and thus indistinguishable, they are shown as 
one prediction in Fig. [TJ The Mif-scheme deviates very slightly from the other schemes 
especially for small tp, this difference is at most of order 0.3%. However in this case the 
uncertainty introduced by the Mif -scheme is an order of magnitude smaller compared to 
the total correction which is of order 3 — 4%. For the sbottom, the corrections are due 
essentially to the QCD/SQCD corrections increasing with tp from 3% to about 4%. This 
correction is by itself small. The correction in the case of the stau mass is even smaller 
by an order of magnitude at least. However, here the MH uncertainty at small tp is 
noticeable at small tp of order 0.1 — 0.2% from the other three tp schemes which agree with 
each other to better than 0.01%. The reason for the (almost) scheme independence is that 
the t^-scheme dependence of the sbottom mass as well as of the stau mass is proportional 
to s\p ~ A/ti which is strongly suppressed for large tp. Our results for the DCPR and 
DR schemes are in excellent agreement with those of Ref. [33]. Concerning the choice of 
the mixing parameter Smj i2 , we observe a small difference between the default choice in 
SloopS given by Eq. (I3.33P and the one given by Eq. (I3.34p . To give an idea, the difference 
is about 0.2% in the sbottom mass correction for both tp = 10 and tp = 50. 
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Figure 1: Heaviest sbottom mass, mg and heaviest stau mass, , at tree-level (solid) 
and at one-loop for the A TT ( and also, DR and DCPR) scheme (dashed) and for the MH 
scheme (dash- dot- dotted) as a function oft/3. The percentage correction is also given. 



5.2 Corrections to the masses of the heaviest neutralinos, m v o 

' -^2,3,4 

We calculated the quantum corrections to the masses of the three neutralinos for the dif- 
ferent schemes of tp implemented in our code and compared our results with Ref. [25] 
which works within the .DCP-R-scheme but otherwise takes the same input parameters, 
namely the chargino masses and the LSP mass. The input chargino/neutralino param- 
eters are m+ = 180 GeV, m + = 350 GeV, m v o = 160 GeV, we choose as in [25] the 
model corresponding to = e x = 1 in Eqs. (14.46 p . (14.471) in order to reconstruct the three 
fundamental parameters Mi, M 2 and /i. The other input parameters given in Ref. [25] 
which enter indirectly in the loop calculation are M^o = 150 GeV, M 3 = 600 GeV and for 
the sfermion sector, M Ll = M- eR = Mq l = M qR = 300 GeV, A e = A d = 900 GeV and 
A u = 100 GeV. 

Fig. [2] shows our findings. The results obtained by using the DCPR- scheme is in complete 
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Figure 2: Neutralino masses at tree-level (solid/black) and at one-loop by using the A TT - 
scheme (dashed/blue), the DR scheme (dotted/light green), the DC PR-scheme (dash- 
dotted/purple) and the MH-scheme (dash- dot- dotted /red) as a function oft/3. 
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agreement with Fig. 2 of the Ref. [25]. The corrections within the A Tr -scheme, DR-scheme 
and DC P R-scheme are very modest. They are largest for m x o where they reach a maxi- 
mum of 5 GeV, which corresponds to a mere 2.5% relative correction. The results between 
the A^-scheme and the DR-scheme are almost indistinguishable, for all value of tp and all 
three masses. DCP-R-scheme is also very close to the latter schemes, a slight deviation oc- 
curs for values of tp in excess of 30. The largest corrections are found with the Mif -scheme 
which deviates considerably from all other schemes when tp is in excess of 30. Therefore 
once again this scheme does not look very suitable. 



5.3 Some decays of the two charginos 

We compute the full electroweak corrections to a few decays of the charginos that were 
considered in Ref. [31] with the help of the code Grace-SUSY at one-loop. One of the main 
differences between our approach and the one adopted in [M] is the definition of tp. In [31] 
tp is closely related to our MH definition. [M] works with renormalised mixing matrices 
apart from the case of sfermions where a shift in the angle defining the diagonalising matrix 
is performed. To conduct the comparison we take set (A) of [M] given in Tabled moreover 
we have m x + = 184.2 GeV, m x + = 421.2 GeV, m x o = 97.75 GeV. We study also the tp 
scheme dependence of the result. As we find an excellent agreement at tree-level, Table [3] 





M A o 






e 




f 


10.00 


424.90 


399.31 


Mr 


184.12 


184.11 


182.19 


Mi 


M 2 


M 3 


Mr 


118.01 


117.99 


111.29 


100.12 


197.52 


610.00 


M 


-398.93 


-452.58 


-444.84 




u 


d 


c 


s 


t 


b 


M Qr, 


565.97 


565.91 


453.05 




546.78 


544.95 


546.84 


544.97 


460.52 


538.13 




-775.58 


-979.08 


-784.72 


-1025.74 


-535.40 


-938.50 



Table 2: Set of supersymmetric parameters defined as set (A) in [3$ . All mass parameters 
are in [GeV]. 

shows the tree-level result for both codes in one column. For one loop results the agreement 
is generally good when we switch to the MH scheme apart from the corrections to the 
X3 4 masses where a difference is noticeable [§. The correction to the X2 ma ss is quite good. 
The corrections to the masses are negligible especially in the DR scheme and A TT scheme. 
In the one-loop corrections to the decays this additional negligible mass correction is not 

§Note that we have found perfect agreement with Ref. [25] as concerns corrections to all Xi (i=2,3,4) 
masses in the DR and DCPR scheme, see Section [5~2l 
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taken into account in a decay such as xt — > W + x® f° r example especially because of the 
large mass difference between xt an d For the decays, the largest discrepancy is for 
xt — > W + Xi and xt ~~ ¥ %xt- However we note that when this discrepancy is largest, the 
correction within our MH scheme deviates drastically from the prediction within the A TT 
and DR schemes. The MH scheme leads, in some decays, to too large corrections. For 
example for xt ~^ W + Xi the MH scheme gives 23% correction whereas the correction 
in DR is only 5%. A similar observation can be made for xt ~^ ^t v r where in the A TT 
scheme the correction is ~ 0% whereas it reaches 24% within our MH scheme. These 
examples also show that for all decays considered in Table [3] the predictions of the A TT 
and DR are within 2% and very often even much better. Once more these examples show 
that the MH scheme is not to be recommended, we suspect strongly that the differences 
we find between Grace-SUSY and SloopS are essentially due to the peculiar choice of the 
scheme based on the heavy neutral CP-even Higgs that greatly amplifies the corrections 
and the differences. 



Decays [GcV] Tree Level Grace SloopS MH SloopS DR SloopS A TT 
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Table 3: Some X12 decays at tree level and at one-loop with three different tp-schemes in 
SloopS compared to Grace-SUSY for set (A) defined in Table® Corrections to the masses 
of X2,3,4 are a ^ so given. 



5.4 e+e -> XiXi 

We now turn to the full 0(a) correction to chargino production at a linear collider. We 
consider the same process as the one computed in [M] within Grace-SUSY, namely e + e~ —>■ 
xtXiil)- We use the same set of parameters Set (A) defined in Tableland study the energy 
dependence of the total cross section. The same cross section has been studied in [351 ESj • 
The QED radiation in view of an event generator has been studied in [37]. Fig. [3] shows 
the cross section of this process computed at tree-level and also at one-loop for different tp- 
schemes. We find excellent agreement with the results of Ref. [31] when specialising to the 
MH-scheme. The A TT , DR and DCPR give corrections within the per-mil level and one 
can hardly distinguish between the three schemes. For this process and with Set (A), the 
MH scheme gives systematically about —1% to —1.5% difference from the other schemes 
which is very small compared to the discrepancies we have noted for some decays of the 
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charginos with the same set of parameters. This suggests that the tp scheme dependence 
is quite small and explains why our results for this process agree very well with those of 
Grace-SUSY. In any case over the whole range of energies the full 0(a) corrections in the 
DR scheme amounts to about —9% for a centre of mass energy y/s = 500 GeV reaching 
a maximum of about —7% at \/s = 700 GeV and dropping to about —11% at s/s = 1300 
GeV. 
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Figure 3: Total cross section of e + e~ — > XiXiil) as a function of y/s at tree-level 
(solid/black) and at one-loop (full order 0(a) ) in the A TT -scheme (dashed/blue), the DR 
scheme (dotted/light green), the DC PR-scheme (dash-dotted/purple) and the MH -scheme 
( dash- dot- dotted /red). The right panel gives the percentage correction. In the left panel con- 
sidering that the A TT -scheme, the DR and the DCPR-scheme are not distinguishable we 
therefore only show the result of the A TT -scheme beside the tree-level and the MH -scheme. 



e + e~ — > f\fi, f 2 f 2 , fif 2 have been calculated in Ref. [38j [39j HO]. In Ref. [38l [39] only 
the electroweak non QED corrections are computed, the QED corrections are dismissed 
on a diagrammatic level by leaving out one-loop Feynman diagrams with virtual photon 
exchange. In Ref. [10] the full 0(a) is performed with a resummation of the leading 
log QED corrections within a structure function approach for the universal initial state 
radiation. We perform here a complete O(a) calculation of these processes and compare 
our results to those of [39] as concerns the electroweak non QED corrections. We therefore 
take scenario 1 of [39] with the following set of parameters: tp = 20, fi = 1000 GeV, 
Mi = 94.92 GeV, M 2 = 200 GeV, M 3 = 669.18 GeV, M A o = 300 GeV, M Ll = M~ eR = 
M 6 = M~ ,- = 400 GeV, At = -500 GeV, Mr = 360 GeV and Mr = 440 GeV. 

WL u R ,d R ' J ' tfl Or 

In [39] the electromagnetic coupling is not taken in the Thomson limit but is fixed from 
a (M§) with a (M§) = 1/127.934. This absorbs large logarithms compared to our on- 
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shell scheme based on a(0) = 1/137.036. In [39] the mixing parameter in the stau sector 
is parameterised through the mixing angle which is renormalised according to Eq. (13.341) . 
For the sake of comparison we will here also switch to this scheme for the sfermion mixing. 




Figure 4: Total cross section of e + e — > 7^(7) as a function of \/s at tree-level 
(solid/black) and at full one-loop in the DC PR scheme (dashed/blue). We also show 
the tree-level improved cross section with a MS (M^) (dash- dotted/red) and the pure weak 
correction in the on-shell scheme as defined in the text (dotted/purple). The full 0(a) rela- 
tive corrections for the three channels with respect to the tree-level cross sections with a(0) 
is shown in the panel on the right (dashed/blue). We also show the weak non QED relative 
correction (dotted/purple) where the improved tree-level cross sections with a MS (M^) is 
used to absorbs large logs from the running of a. This correction should be contrasted to 
the one obtained in Ref. [3^. In order not to crowd the figure the channels are not labeled. 
They can be easily identified as they have different thresholds. 



In addition to the tree level cross section calculated with a(0) = 1/137.036 and the com- 
plete 0(a) one-loop correction, we compute the improved tree-level cross section with 
a (M§) = 1/127.934. Our evaluation of the weak non QED correction is different from 
the one in [39]. In our case the weak correction is obtained by subtracting the leading 
QED corrections. The initial state radiation factor including the virtual photon correction 
and the soft bremsstrahlung photon below the cut-off energy k c is universal and known, 
see for example |41j . 

%™ = v - 1} ln I + I Le + T ~ ' Le = Hs/ml) ■ (5 ' 1} 

where m e is the electron mass and Eb the beam energy s = 4Ef. To subtract not only the 
initial but also the final state radiation and the final-initial interference QED effect, we 
take the result of the virtual one-loop correction and the soft radiation factor obtained by 
the code and subtract the following 

a wcak (^) = a M+so& (VS, k c ) - -A(Vs) hx(^)- ^a trcc (^) ln ( V (5.2) 

71 \ \/s I 27r V mi I 



36 



The last term in Eq. (15. 2p stems from the collinear singularity due to initial state radiation 
and we neglect non log terms, the latter that arise from initial radiation are negligible of 
order 0.3% relative correction. The term A(a/s) is extracted numerically based on the fact 
that the weak non QED correction is independent of the cut-off k c . We take two small 
enough cut-off k Cl , k C2 to extract A(y/s), 

-.virtual+soft/ l~Z t, \ ^virtual+soft/ fZ h, \ 

-A(y/s) = Ws,k Cl )-a . (5.3) 

We have checked that a weak (y/s) defined this way is independent of the cut-off k c by taking 
other values of k r . Such a definition of the weak correction has been introduced in [421 . 



Our tree-level results for the improved tree-level with a = a MS (M^) reproduces the corre- 
sponding cross section in Ref. [39] perfectly. 

To help compare our results with those Ref. |39j . the right panel of Fig. H] shows also 
the relative weak non QED correction with a M (Mf) as input rather than a(0), hence 
subtracting large logs from the running of a. Our predictions for the weak correction de- 
fined this way are within 1% of those in [39] within the DCPR scheme used in [39]. We 
have traced this small difference to the different ways the weak correction is defined from 
the subtraction of the QED corrections. The energy dependence of the weak corrections 
matches perfectly. 

We can now comment on the tp scheme dependence and the sfermion mixing renormalisa- 
tion scheme. The corrections induced by the different tp schemes are very small. Even the 
MH scheme departs by not more than 0.3% from the DR. The other schemes, DCPR 
and A TT , agree within better than 0.01% with DR. The difference in the choice of the 
sfermion mixing parameter Srnj i2 is even more negligible here. For example, for a centre of 
mass energy yfs = 1000 GeV, the one-loop correction to the process e + e~ — > f\fi differs 
only about 0.003% when we switch from the default definition in SloopS Eq. (13.331) to the 
one that has been usually used Eq. (13.341) . 

These calculations show that not only it is important to take into account the QED cor- 
rections but also that the pure electroweak corrections are certainly not negligible, for 
example even after absorbing the effect due to the running of a, the weak corrections for 
f\fi production is about —15%. 



6 Conclusions 

We have presented in detail a complete renormalisation of the sfermion sector as well as of 
the chargino/neutralino sector of the MSSM in the case of CP conservation. We critically 
analysed the renormalisation of the mixing parameter in the sfermion sector and discussed 
different ways to define it in a consistent manner. This paper is a sequel to our study in 
Ref. [1] and completes the presentation of all the ingredients that are built into our au- 
tomatised code for one-loop calculations in the MSSM, SloopS. Although other approaches 
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to renormalising the MSSM have been worked out, we believe that our approach treats 
all the sectors consistently within the same general on-shell framework in particular about 
the treatment of mixing and how one deals with the rotation and diagonalising matrices. 
Moreover our code permits powerful gauge checks with the help of the non-linear gauge fix- 
ing condition and allows to easily switch between different renormalisation schemes. Some 
very powerful and extensive tests have been conducted on the code as concerns ultraviolet 
finiteness and gauge parameter independence on an almost exhaustive list of 2 — > 2 pro- 
cesses, see [31]. In the present paper we choose to concentrate on a few key observables 
in the sfermion and chargino/neutralino sector and compared our results with some that 
are found in the literature while at the same time studying the impact of different renor- 
malisation schemes. We have calculated one-loop corrections to sfermion masses and also 
neutralino masses. We have also derived some chargino decay widths and presented a cal- 
culation of the production of charginos and sleptons at e + e~ colliders. We find the genuine 
electroweak corrections in these cross sections to be rather important and should therefore 
be taken into account. Having at our disposal a code that allows the one-loop calculation 
for any process in the MSSM, it is now possible to envisage revisiting analyses for the 
extraction of the fundamental supersymmetric parameters from precision measurements 
at the colliders and use them in turn for a precision calculation of the relic density for 
example. Finally, let us mention that other renormalisation schemes, with different choices 
of the input parameters from the one described in this paper, for the chargino/neutralino 
sector are already implemented in the code and would be part of a forthcoming study. 
Although in the many examples we have shown here the QCD corrections are calculated, 
a complete treatment of the gluon/gluino sector within an automated code such as SloopS 
and in particular how to easily implement within the code a regulator for the infrared 
singularity is work in progress. 



Acknowledgments 

We would first like to thank Andrei Semenov whose help was invaluable in the first stages 
of the project. We also owe much to our friends of the Minami-Tateya group and the devel- 
opers of the Grace-SUSY code, in particular we learned much from Masaaki Kuroda. We 
benefited a great deal from discussions with Guillaume Chalons, Sun Hao, Karol Kovarik 
and Peter Zerwas. This work is supported in part by GDRI-ACPP of the CNRS (France). 
This work is part of the French ANR project, ToolsDMColl. This work is also supported 
in part by the European Community's Marie-Curie Research Training Network under con- 
tract MRTN-CT-2006-035505 "Tools and Precision Calculations for Physics Discoveries 
at Colliders", the DFG SFB/TR9 "Computational Particle Physics", and the Helmholtz 
Alliance "Physics at the Terascale". 



38 



References 



[1] N. Baro, F. Boudjema, A. Semenov, Phys. Rev. D78 (2008) 115003. larXiv:0807.4668l 
[hep-ph] . 

[2] A. Semenov, |hep-ph/9608488] 

A. Semenov, Nucl. Inst. Meth. and Inst. A393 (1997) 293; 
A. Semenov, Comp. Phys. Commun. 115 (1998) 124; 



A. Semenov, 
A. Semenov, 



hep-ph/0208011; 



arXiv:0805.05 55 [hep-ph]. 



[3] J. Kiiblbeck, M. Bohm, A. Denner, Comp. Phys. Commun. 60 (1990) 165; 
H. Eck, J. Kiiblbeck, Guide to FeynArts 1.0, Wiirzburg, 1991; 
H. Eck, Guide to FeynArts 2.0, Wiirzburg, 1995; 
T. Hahn, Comp. Phys. Commun. 140 (2001) 418, hep-ph/ 00 12260. 



[4] T. Hahn, M. Perez- Victoria, Comp. Phys. Commun. 118 (1999) 153, hep-ph/9807565 



T. Hahn, hep-ph/0406288 hep-ph/0506201 



[5] T. Hahn, LoopTools, http://www.feynarts.de/looptools/. 

[6] N. Baro, F. Boudjema, A. Semenov, Phys. Lett. B660 (2008) 550, larXiv:0710.182T1 
[hep-ph]. 



[7] F. Boudjema, A. Semenov, D. Temes, Phys. Rev. D72 (2005) 055024, hep-ph/0507127 



[8] J. Fujimoto et al, Comput. Phys. Commun. 153 (2003) 106, hep-ph/0208036 



[9] G. Belanger, F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato, Y. Shimizu, 
Phys. Rep. 430 (2006) 117, |hep-ph/0308080[ 



[10] F. Boudjema, E. Chopin, Z. Phys. C73 (1996) 85, hep-ph/9507396 



[11] A. Freitas, D. Stdckinger, Phys. Rev. D66 (2002) 095014, hep-ph/0205281 



[12] J. A. Coarasa, D. Garcia, J. Guasch, R.A. Jimenez, J. Sola, Eur. Phys. J. C2 (1998) 



373, hep-ph/9607485 



J. A. Coarasa, D. Garcia, J. Guasch, R.A. Jimenez, J. Sola, Phys. Lett. B425 (1998 
2 

Guasch, J 



329, |hep-ph/9711472 
J.A 



Coarasa, J 
|hep-ph/9808278| 

[13] P.H. Chankowski, S. Pokorski, J. Rosiek 



Sola, W. Hollik, Phys. Lett. B442 (1998) 326, 
Nucl. Phys. B423 (1994) 437, 



hep-ph/9303309, 



[14] A. Denner, Fortsch. Phys. 41 (1993) 307. larXiv:0709.1075l [hep-ph]. 



[15] A. Arhrib, R. Benbrik, Phys. Rev. D71 (2005) 095001, hep-ph/0412349 



[16] A. Bartl, H. Eberl, K. Hida ka, S. Kraml, W. Majerotto, W. Porod, Y. Yamada, Phys. 
Lett. B419, 243 (1998), [hep^ph/9710286[ 

[17] R.G. Stuart, Phys. Lett. B262 (1991) 113. 



39 



[18] J.R. Espinosa, I. Navarro, Phys. Rev. D66 (2002) 016004, [hip^ph/0109126 



[19] J.R. Espinosa, Y. Yamada, Phys. Rev. D67 (2003) 036003, hep-ph/0207351 



[20] J. Guasch, J. Sola, W. Hollik, Phys. Lett. B437 (1998) 88, hep-ph/9802329 



H. Eberl, S. Kraml, W. Majerotto, JEEP 9905 (1999) 016, hep-ph/9903413 



[21] S. Kraml, PhD dissertation, hep-ph/9903257 



[22] H. Eberl, A. Bartl, W. Majerotto, Nucl. Phys. B472 (1996) 481, |hep-ph/9603206t 
S. Kraml, H. Eberl, A. Bartl, W. Majerotto, W. Porod, Phys. Lett. B386, 175 (1996), 



hep-ph/9605412; 



A. Bartl, H. Eberl, K. Hidaka, S. Kraml, W. Majerotto, W. Porod, Y. Yamada, Phys. 
Lett. B419, 243 (1998), |hep^ph/9710286[ 

A. Bartl, H. Eberl, K. Hidaka, S. Kraml, W. Majerotto, W. Porod, Y. Yamada, Phys. 



Rev. D59, 115007 (1999), hep-ph/9806299 



[23] A. Djouadi, W. Hollik, C. Junger, Phys. Rev. D55, 6975 (1997), [hep^ph/9609419 

[24] J. Guasch, J. Sola, W. Hollik, Phys. Lett. B437 (1998) 88. 

[25] T. Fritzsche, W. Hollik, Eur. Phys. J. C24 (2002) 619, |hep-ph/0203159| 



[26] M. Drees, W. Hollik, Q. Xu, JEEP 02 (2007) 032, hep-ph/0610267 



[27] J. Guasch, W. Hollik, J. Sola, JEEP 0210 (2002) 040, hep-ph/0207364 



[28] J.L. Kneur, G. Moultaka, Phys. Rev. D59 (1999) 015005, hep-ph/9807336 



[29] S. Y. Choi, A. Djouadi, M. Guchait, J. Kalinowski, H. S. Song, P. M. Zerwas, Eur. 
Phys. J. C14 (2000) 535, |hep-ph/0002033| 

[30] SY. Choi, J. Kalinowski, G.A. Moortg at-Pick, P.M. Zerw as, Eur. Phys. J. C22 (2001) 
563, Addendum-ibid. C23 (2002) 769, |hep-ph/0108117| 

[31] N. Baro, PhD thesis, Renormalisation and predictions at one-loop in super symmetry, 
applications to dark matter and collider physics, 



http : //tel . archives-ouvertes . f r/tel-00329722/f r/ 



[32] [CompHE P Collab oration], E. Boos et al, Nucl. Instrum. Meth. A534 (2004) 250, 
|hep-ph/0403113[ 

A. Pukhov et al, "CompHEP user's manual, v3.3", Preprint INP MSU 98-41/542 (1998) 
|hep-ph/9908288t 

http : //theory . sinp . msu . ru/ comphep/. 



[33] W. Hollik, H. Rzehak, Eur. Phys. J. C32 (2003) 127, hep-ph/0305328 



[34] J. Fujimoto, T. Ishikawa, Y. Kurihara, M. Jimbo, T. Kon, M. Kuroda, Phys. Rev. 
D75 (2007) 1 13002, [hep^ph/0701200 



[35] T. Fritzsche, W. Hollik, Nucl. Phys. Proc. Suppl. 135 (2004) 102, |hep-ph/0407095[ 



[36] W. Oiler, H. Eberl, W. Majerotto, Phys. Rev. D71 (2005) 115002, hep-ph/0504109 



[37] W. Kilian, J. Reuter, T. Robens, Eur. Phys. J. C48 (2006) 389, hep-ph/0607127 



40 



[38] A. Arhrib, W. Hollik, JEEP 0404 (2004) 073, |hep-ph/03lil~49 



K. Kovarik, C. Weber, H. Eberl, W. Majerotto, Phys. Lett. B591 (2004) 242, 



hep-ph/0401092 



[40] K. Kovarik, C. Weber, H. Eberl, W. Majerotto, Phys. Rev. D72 (2005) 053010, 



hep-ph/050602lj 



[41] G. Belanger, F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato, Y. Shimizu, 



Phys. Lett. B559 (2003) 252, |hep-ph/02 12261 



[42] N. Bouayed and F. Boudjema, Phys. Rev. D77 (2008) 013004. [irXiv:0709.4388l [hep- 
ph]. 



41 



